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ABSTRACT 

We present HSiM: a dedicated pipeline for simulating observations with the HARMONI 
integral field spectrograph on the European Extremely Large Telescope. HSIM takes 
high spectral and spatial resolution input data-cubes, encoding physical descriptions 
of astrophysical sources, and generates mock observed data-cubes. The simulations 
incorporate detailed models of the sky, telescope and instrument to produce realistic 
mock data. Eurther, we employ a new method of incorporating the strongly wave¬ 
length dependent adaptive optics point spread functions. HSIM provides a step beyond 
traditional exposure time calculators and allows us to both predict the feasibility of 
a given observing programme with HARMONI, as well as perform instrument design 
trade-offs. In this paper we concentrate on quantitative measures of the feasibility of 
planned observations. We give a detailed description of HSIM and present two studies: 
estimates of point source sensitivities along with simulations of star-forming emission¬ 
line galaxies at 2 : ^2-3. We show that HARMONI will provide exquisite resolved 
spectroscopy of these objects on sub-kpc scales, probing and deriving properties of 
individual star forming regions. 

Key words: instrumentation: spectrographs - instrumentation: adaptive optics - 
galaxies: high-redshift - galaxies: kinematics and dynamics 


1 INTRODUCTION 

As the current suite of 8 m class telescopes are being pushed 
to their limits, it is becoming more important to prepare for 
and define the scientific programmes to be undertaken by 
the next generation of extremely large telescopes. The Euro¬ 
pean Extremely Large Telescope (E-ELT) is a 39 m diameter 
telescope being constructed by the European Southern Ob¬ 
servatory (ESO) for operations commencing within the next 
decade. Upon completion, this facility will provide the angu¬ 
lar resolution and light gathering power to revolutionise our 
current understanding in many areas of astrophysics. How¬ 
ever, as telescope and instrumentation projects enter such 
ambitious scales of size and complexity of design, it is be¬ 
coming imperative to accurately quantify performance prior 
to construction. 

Consequently, several recent integral field spectrograph 
(IFS) instrumentation projects have developed detailed 


* E-mail: simon.zieleniewski@physics.ox.ac.uk 


instrument simulation models, including KMOS (|Lorente| 
et al.||20Q8|), MUSE (jJarno et al ||2QQ8t , JWST/NIRSpec 
(Piqueras et al. 2010 Dorner et al.|2011 ), and EAGLE/ELT- 


MOS ( Puech et al.|2008 2010 Puech, Yang Flores|2010| . 

The most rigorous simulation method involves using the full 
optical design of the instrument to propagate photons to 
the detector plane, before using data reduction software to 
reduce the simulations akin to real data. While this full end- 
to-end technique is possible for currently operational instru¬ 
ments and those at an advanced stage of design, it is a level 
above what is required at an early stage of development. 


As part of the concept design for an ELT multi-object 
spectrograph (MOS), [Puech, Yang Flores (2010) devel¬ 
oped an IFS instrument simulator that processed input 
data-cubes, generating mock observation cubes assuming a 
perfect data reduction process. The simulator encoded all 
the sky, telescope and adaptive optics (AO; through the 
point spread function) and instrument parameters, return¬ 
ing mock observations containing source and background 
flux as well as noise contributions. They were able to anal- 
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yse the feasibility of a key MOS program to observe star 
forming galaxies at z ~ 6. 

The High Angular Resolution Monolithic Optical and 
Near-infrared Integral field spectrograph (HARMONI) is a 
proposed instrument selected as one of two first-light in¬ 
struments by ESO ( Thatte et al.|20fQ ). HARMONI is being 
designed as a single-field, visible and near infra-red (NIR) 
IFS, and will provide a range of spatial pixel (spaxel) scales 
and spectral resolving powers, which permit the user to op¬ 
timally configure the instrument for a wide range of sci¬ 
ence programs; from ultra-sensitive to diffraction limited, 
spatially resolved, physical (morphology), chemical (abun¬ 
dances and line ratios) and kinematic (line-of-sight veloc¬ 
ities) studies of astrophysical sources. HARMONI will be 
compatible with two modes of adaptive optics (AO) allowing 
it to tackle a broad range of problems across astrophysics, 
including: a) the physics of mass assembly of galaxies at high 
redshifts, b) resolved studies of stellar populations in distant 
galaxies, and c) detecting and weighing intermediate mass 
black holes in nearby galaxies or globular clusters. The full 
design specihcations are listed in table 1. The instrument 
will commence its preliminary design phase in 2015 and so 
here we undertake detailed simulations of the performance 
of HARMONI. 


In this paper we present our development of the instru¬ 
ment simulation pipeline HSIM for the HARMONI instru¬ 
ment that will be used to quantify the performance. HSIM 
improves on the method developed by Puech et al. in sev¬ 
eral key areas, including spectral line spread convolutions, 
and a wavelength dependent telescope (plus adaptive optics) 
point spread function model. We detail the components of 
HSIM and show the importance of the key stages and their 
effects on observations. We then showcase the pipeline with 
two studies: calculations of point source sensitivities for the 
available observing modes, and a case study of measuring 
kinematics of star forming galaxies at z ~ 2-3. 

The paper is organised as follows: Sections and de¬ 
tail the motivation, goals and methodology of the simula¬ 
tions; in Section [4]we present an overview of the pipeline 
stages; in Section^ we present our point spread function 
parameterisation and demonstrate its importance; Section 
shows our verification crosschecks against existing software; 
in Section we present predicted point source sensitivities; 
Section [^presents our high redshift galaxy simulations; and 
we finally conclude in Section 

Throughout this paper we adopt a flat ACDM cosmol¬ 
ogy with Hq — 70kms“^ Mpc“^, Qrn = 0.3 and Qa = 0.7. 
We use AB magnitudes throughout unless otherwise stated. 


2 MOTIVATION AND GOALS 

Developing an instrument simulation pipeline serves several 
important purposes: 

a) it gives a quantitative understanding of the feasibility 
of observing programs, characterising both the capabilities 
and limitations of an instrument, 

b) it allows for performance trade-offs between differing 
instrument designs or configurations. 

The first point is important since it will allow observing 
time on the E-ELT to be more efficiently utilised. Quanti¬ 


fying performance usually consists of determining the ex¬ 
posure time required to achieve a particular signal-to-noise 
ratio (SNR) for a given target. However, an instrument simu¬ 
lation pipeline allows us to go one step beyond exposure time 
calculators and ascertain the precision with which we can de¬ 
rive a number of key physical parameters for particular sci¬ 
ence cases. The input data-cubes can be built to encode as 
many physical characteristics as desired. The mock output 
data represents that of real observations and so it is possible 
to use identical analysis methods as used on real data (or 
permits users to develop and verify analysis methods prior 
to taking the real data). This allows the user to specifically 
determine the feasibility of their science goal (e.g. measuring 
kinematics, black hole masses, stellar chemical abundances). 
It also provides a direct comparison between the scientihc 
capabilities of current generation instrumentation and those 
of the E-ELT era. 

Designing a self-contained and modular pipeline means 
that all simulations are kept consistent within the current 
design of the instrument. The pipeline can easily be modified 
to reflect any changes in the instrument design (e.g. changing 
spaxel scales, spectral resolving power etc), thus allowing ex¬ 
ploration of the second point above. It will also be possible to 
explore instrument effects such as imperfect sky subtraction 
and spectrograph throughput variations, and give quanti¬ 
tive constraints on such calibration issues. In this paper we 
concentrate on the accuracy with which physical parameters 
can be measured for a specific science case. 


3 SIMULATION METHODOLOGY 

The procedure for performing simulations can be broken 
down into three parts: 

a) An input data-cube encoding the physical characteristics 
of the object. This can include properties such as kine¬ 
matics, absorption/emission lines, dynamics, morphology, 
chemical abundances etc; 

b) A simulation pipeline that takes the input and adds 
all the first-order sky, telescope, instrument, and detector 
effects, as well as random and systematic noise, creating 
an output mock-observed data-cube. The user is able to 
choose suitable observing parameters to be able to explore 
the instrument and telescope modes; 

c) The analysis of the output mock-observed data to enable 
an understanding of what information can be extracted 
from the observation, i.e. how well can the properties of the 
input data-cube be recovered in the output cube. 

Within this framework there are two areas of parame¬ 
ter space to be explored: the physical characteristics of the 
objects (e.g. morphology, mass, redshift, magnitude, line 
strengths) as well as the observational characteristics of the 
telescope, and instrument (e.g. spatial scale, exposure time, 
resolving power, emissivity, AO correction). Consequently 
the combined parameter space is vast, and a full exploration 
is simply not possible. For the purposes of this paper, we 
have firstly pursued point source sensitivity calculations for 
all current instrument configurations, which gives a general 
overview of the performance of the instrument design. We 
then present simulations of a representative case oi z ^ 2-3 
emission-line galaxies. 
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Table 1. Current design parameters for HARMONI. 

Wavelength range 0.47 - 2.45 p,m 

Spatial scales (FoV) 4 mas (0.86" x 0.61"), 10 mas (2.14" x 1.52") 

20 mas (4.28" x 3.04"), 60 x 30 mas (9.12" x 6.42") 

Spectral Resolution (& waveband) 3500 {V + R; Iz J; H K) 

7500 (/z; J; H; K) 
20000 (half of each NIR band) 

120 K cryostat, 40 K detectors 

Teledyne HgCdTe 4k x 4k for NIR 
4k X 4k CCD for visible 

Throughput Target ^ 30% average 

AO modes SCAO, LTAO, seeing-limited 


Temperature 

Detectors 


Input data-cube 

Spatial structure 

Velocity scale 
Emission/absorption lines 

Spectral continuum 


Observing parameters 

DIT, NDIT, grating, spaxel scale, 

AO mode, seeing, Z^, 


Spectral LSF convolution 



Atmospheric differential refraction 


Spatial PSF convolution 
Spaxel rebinning 

Sky, telescope, instrument, detector 
background & throughput models 


Poisson and systematic noise 



Figure 1. Flowchart of the simulation process. Input data-cubes 
provide all the physical details of the object. The cube is then ‘ob¬ 
served’ for given instrument, telescope and site parameters. The 
observing process adds all first-order telescope and instrument ef¬ 
fects, as well as random and systematic noise. The output cube 
represents a perfectly reduced data-cube, and can be analysed 
exactly like real data. 


4.1 Line spread function convolution 

The spectral resolution of the input cube is degraded to the 
output resolution depending on the chosen grating. This is 
achieved by convolving the spectral dimension with a Gaus¬ 
sian line-spread function (LSF). We use a Gaussian of con¬ 
stant full width half maximum (FHWM) as a reasonable ap¬ 
proximation for the resolution of a slit-width limited grating 
spectrometer. We ensure the output resolution of the cube 
is consistent by generating a LSF of width AAconv given by 

AAconv = y(AAout)2 - (AAi„)2, (1) 

where AAin is the spectral resolution of the input cube speci¬ 
fied in the fits header (by the user), and AAout is the output 
resolution given by the chosen grating. The spectral dimen¬ 
sion of the cube is then sampled (with a minimum of two 
pixels per FWHM to approximate the Nyquist limit). If the 
input spectral resolution is coarser than the grating reso¬ 
lution (AAin > AAout) then the user can choose to either 
convolve with a LSF of width AAout or perform no spectral 
convolution. The former option is useful in situations where 
the input data-cube contains a spectrum where the resolu¬ 
tion is limited by e.g. stellar velocity dispersion in galaxies, 
and the spectrum still needs to incorporate the instrumental 
effect. The option to ignore the spectral convolution is al¬ 
ways available if it is not required (e.g. if the input data-cube 
has been pre-formatted spectrally). 


4 SIMULATION PIPELINE OVERVIEW 


We first presented the simulation pipeline in |Zieleniewski| 
et ah (2014) and here we provide an overview of each stage, 


highlighting several novel features. The overall process is 
presented in Fig. HSIM is written in python and makes 


use of the astropy package (Astropy Gollaboration et al. 


2013), as well as the numpy and SCIPY packages. It han¬ 
dles input and output data in fits format, as this is the 
usual file format for 3D spectrographs. The input data-cube 
is uploaded as a 3D array of two spatial and one spectral 
dimension. The fits headers contain all the relevant infor¬ 
mation about the object (e.g. spatial and spectral sampling, 
flux units, spectral resolution). The user chooses suitable 
observing parameters and the data-cube is then processed 
through several steps: 


4.2 Atmospheric differential refraction 


The effect of atmospheric differential refraction (ADR) is 
added according to the equations of [Schubert Walter-] 


scheid (2000) and Roe (2002). The angle in radians between 


the true zenith distance Zd and the apparent zenith distance 
Za is approximated by 


R — Zd — Za — 


- 1 
2n‘^ 


tan Zd, 


( 2 ) 


where the refractive index n is a function of wavelength, 
pressure, temperature and humidity. Thus the angular 
change Ri — R 2 , between two wavelengths is given by: 

where n is the refractive index at a given wavelength and 
Zd the zenith distance. This effect is prominent at both 
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visible wavelengths, where the refractive index varies the 
most strongly with wavelength, and at the smallest spatial 
scales. Each wavelength channel (spatial image at a given 
wavelength in the 3D cube) is shifted relative to an optimal 
wavelength, which is calculated to give equal shift on either 
side. The shift is made along the longest spatial dimension 
of the output data-cube (this assumes the user aligns the 
longest spatial axis of the instrument field of view along the 
parallactic angle). After the data-cube is rebinned to the 
chosen spaxel scale, the ADR is ‘corrected’ by shifting the 
wavelength channels back again to emulate the correction 
achieved by a data reduction pipeline. This stage reduces 
the common field of view of the data-cube. The second or¬ 
der effect whereby the axis of ADR moves throughout an 
exposure causing the object to blur in the data-cube, is not 
included. 


4.3 Point spread function convolution and spatial 
rebinning 

Each wavelength channel is convolved with an AO spatial 
PSE. We have developed a novel method for parameteris- 


ing AO PSEs as detailed in Zieleniewski & Thatte (2013) 


allowing us to generate PSEs at any wavelength within the 
HARMONI range. This is described in more detail in Sec¬ 
tion The pipeline generates PSEs at Imas sampling and 
performs the PSE convolution at a scale of 1/lOth of the 
output spatial scale in order to minimise convolution effects 
due to the finite sampling of the input cube and PSE. If the 
input data-cube is coarser than 1/lOth of the output scale 
then it is interpolated up to 1/lOth, ensuring flux conserva¬ 
tion. 

The data-cube is rebinned to the chosen output spatial 
scale ensuring flux conservation. 


4.4 Background and throughput 


Total background and throughput cubes are generated in¬ 
corporating the effects of sky, telescope, instrument and de¬ 
tector. We use the ESO Skycalc sky model described by 


Noll et al. (2012) and Jones et al. (2013) for both emission 


and transmission. We model the telescope as a grey-body 
whereby a thermal black-body curve for the given site tem¬ 
perature is multiplied by a constant emissivity. 


4.5 Noise 


Poisson noise from the object, sky and telescope background, 
and detector dark current is added along with detector read¬ 
out noise. We use separate statistics to represent the use 
of CCDs for visible wavelengths and NIR arrays for longer 
wavelengths. The NIR detector statistics are modelled on 
the KMOS HAWAII-2RG detectors ( Einger et al.|20d^ and 
the CCD statistics based on the MUSE calibration data (D. 
Ives, priv. comm.). 


4.6 Outputs 

The final outputs of HSIM are: 

1. A mock observed cube: the main product of the simulation 
containing flux from the source, background and detector 



Figure 2. Strehl ratios of parameterised PSPs used in HSiM as a 
function of wavelength for SCAO (solid lines) and LTAO (dashed 
lines). PSFs are generated at a seeing of 0.67" and sampling of 
1 mas. Both the smooth variation of AO performance and the 
improved performance at longer wavelengths are evident. 


along with all associated noise for each (x, y. A) pixel of the 
cube; 

2. A background cube containing all background flux; 

3. A SNR cube giving the signal-to-noise ratio for each pixel 
calculated as 

_ _ Q(x, y, A)Texp y/Wxp _ 

>/0(x, y, A)Texp + B{x, y, A)Texp + D • Texp + cr^ ’ 

( 4 ) 

where 0{x,y,X) is the object counts per second for a pixel 
in the 3D cube, B(x,y,X) is the total background counts 
per second, D is the total dark current per second, ctr is 
the read-out noise, Texp is the exposure time in seconds and 
Wxp is the total number of exposures. 

There are also options to return a noiseless object cube con¬ 
taining only source flux, return a transmission cube, and 
perform automatic perfect sky subtraction (observed cube - 
background cube). 


5 EFFECTS OF THE AO POINT SPREAD 
FUNCTION 

The effects of the telescope and AO system on observations 
are encoded in the spatial PSE. The E-ELT PSE is predicted 
to be a strong function of wavelength as well as other pa¬ 
rameters including seeing, off-axis distance from reference 
AO star, and guide star magnitude. Its complex form means 
it needs to be incorporated carefully into the simulations 
not simply modelled as a simple Gaussian or Moffat func¬ 
tion. We have incorporated a continuously varying AO PSE 
as a function of wavelength and seeing, for both LTAO and 
SCAO, in our simulation pipeline, using the eltpsffit program 
developed by J. Lisk^ This allows the user to fit a ID aver¬ 
age radial profile of a PSE with a set of analytical functions. 
We used a set of simulated long-exposure LTAO and SCAO 


^ https://www.eso.org/sci/facilities/eelt/science/drm/ 
tech_data/ao/psf.fitting/ - last accessed 01-06-15 
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Figure 3. Integrated continuum maps from an HSIM simulation of a ULIRG observed with LTAO. Left panel shows the input map at a 
sampling of 4 mas after scaling an HST WCS image of a local galaxy (50 mas sampling) to the correct spatial extent at 2 ; = 2. Centre 
panel shows a 20 mas continuum map at 1.0 p,m in the /z-band. Right panel shows a continuum map at 2.2 p,m in the RT-band. The 
observed maps are shown at SNR = 10. We assume identical intensity and spatial distribution for both continuum maps. The improved 
PSF performance at longer wavelengths is clearly evident. 


PSFs (T. Fusco & N. Schwartz, priv. comm.), covering the 
HARMONI wavelength range, and fitted each of these with 
a combination of an obscured Airy function, Moffat function 
and Lorentz function ( Zieleniewski Thatte|2013 ). Interpo¬ 
lating the parameters of each analytical function with wave¬ 
length, we then obtain the PSF at each wavelength. We have 
so far only used on-axis PSFs and we assume the PSF does 
not vary spatially within the field of view. Fig. shows the 
measured Strehl ratios of our parameterised PSFs for both 
LTAO and SC AO. The smooth variation of AO performance 
with wavelength in our parameterised PSFs is clear. 

The importance of the wavelength variation is worth 
emphasising. Whilst the diffraction limit of a telescope in¬ 
creases from visible to NIR wavelengths, AO performance 
vastly improves at longer wavelengths because the atmo¬ 
spheric turbulence becomes relatively less disruptive. This 
large improvement, as seen in Fig.|^ more than compensates 
for the increasing diffraction limit at longer wavelengths. We 
demonstrate this qualitatively in Fig. which shows two 
continuum maps of a simulated HARMONI observation of 
an ultra luminous infra-red galaxy (ULIRG) at a redshift of 
z — 2. The input HST ACS image at a resolution of 50 mas 
(IRAS06076-2139: from Armus et al.||2009| has been scaled 
spatially to 4 mas at 2 : = 2. The left map shows the object as 
observed in the /z-band at 1.0 /xm. The right map shows the 
object as observed in the iF-band at 2.2 /xm. For easier com¬ 
parison we assume equal flux at both wavelengths and we 
ignored the increased thermal background in the iF-band. 
It is clear from the two maps that the spatial resolution in 
the iF-band is much improved over the /z-band, providing 
exquisite detail of the underlying star-forming morphology 
of the ULIRG. 


The AO PSFs that we used for the parameterisation 
also include a small amount of intrinsic jitter or blur, whicht 
estimates the effect of wind-shake jitter on the telescope 
structure. This corresponds to 2 and 3 mas rms for the SGAO 
and LTAO PSFs respectively. We include the ability to add 
additional amounts of PSF blur to approximate reduced AO 
performance. We model this by convolving the AO PSF with 
a Gaussian of chosen rms width, so any additional blurring 
adds in quadrature. 


6 VERIFICATION RUN 


To test the results of HSIM we compare to two separate tools. 


Firstly, we compared to the pipeline of Puech et ah (2010), 


which has itself been verified using real SINFONI observa¬ 
tions from Genzel et al. ( |2006 ). As part of the E-ELT design 
reference mission (DRM) science case GIO: The physics and 
mass assembly of galaxies out to z ^ 6, |Puech et al.| under- 
took simulations of which some examples are publicly avail¬ 
ably We took the example data-cube UGG5253, a rotating 
disc galaxy with prominent Ha emission. 

We converted the cube for compatibility with our 
pipeline using the description from the DRM. We then ran 
simulations using identical parameters as used by |Puech| 
|et al.| including the identical spatial PSE. The only differ¬ 
ence between the two runs is that Puech et al. do not spec¬ 
ify the sky transmission spectrum they use, so we manually 
convolved one from the E-ELT DRM to match the spectral 
resolution of the data-cube. 

The outputs of the comparison are shown in Eig.^ Our 
pipeline computes the total signal of iVexp combined expo¬ 
sures, so we have divided our output cube by iVexp to get a 
mean value for each pixel (cf. Puech et al. who determine 
the median value of each pixel from Vgxp cubes). It is clear 
from both the signal and SNR plots that we are computing 
consistent values within the uncertainty of different trans¬ 
mission functions, and can be confident that our pipeline is 
working correctly. 

Eor the second crosscheck we compared our pipeline 
with the ESO SINEONI exposure time calculator (ETG). We 
ran the ETG for an extended object with an AOV spectrum 
at a iF-band Vega magnitude of 14, observed with no adap¬ 
tive optics for two hours (8 x 900 s). The ETG returns the 
input spectrum which we used to create an input cube. We 
simulated an identical observation with HSIM, using a Gaus¬ 
sian PSE and the relevant VLT telescope parameters. We 
were careful to calculate our object flux over two spatial pix¬ 
els as done by the ETG. We note that we approximated the 


^ https://www.eso.org/sci/facilities/eelt/science/drm/ 
CIO/ - last accessed 01-06-15 
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Figure 4. Left: comparison of the signal extracted from a 10 spaxel radius aperture centred on the galaxy between Puech et al. simulation 
and HSIM. Right: comparison of the median SNR over every spaxel in each wavelength channel. The simulations are consistent to within 
the noise realisations. 


output resolution of our spectrum using the central wave¬ 
length value of the spectrum divided by the resolving power 
of the SINFONI grating, as it is unclear exactly how the 
LSF convolution is performed in the ETC. We show the re¬ 
sults of a comparison in Fig. The two main panels show 
the total SNR and object signal within two spaxels. The 
two smaller panels show the respective residuals after me¬ 
dian smoothing by 10 spectral pixels to account for potential 
differences in LSF convolution. Our results show extremely 
good agreement with the ETC with sub 5 per cent residu¬ 
als over the majority of the spectrum (the largest residuals 
appearing at the blue end of the spectrum where the flux 
values are very small and edge effects occur). This gives us 
further confidence that our pipeline is working correctly and 
is producing realistic signal and noise values. 


7 SENSITIVITY PREDICTIONS 


In this section we present instrument point source sensitiv¬ 
ities for several operating configurations of HARMONI as 
computed using HSIM. We compute the limiting AB mag¬ 
nitude Mab to achieve SNR = 5 per spectral resolution 
element (2 pixels) after five hours on source (Texp = 900 s). 
We assume sky subtraction with no penalty to on-source ex¬ 
posure time, using ‘nodding-on’ or similar techniques. To do 
this we calculate the SNR for a given Mab after 5 hours on 
source, then compute the multiplicative factor / required 
to give a Mab that yields SNR = 5. Adding this factor 
(and extra factors for sky subtraction and summing over 
two spectral pixels) equation becomes 


SNR = 


f ‘ O ‘ Texp Aexp 

\/(/ • O + 2B + 2D) • Texp + 2(7^ 


( 5 ) 


and this is solved for / using the standard quadratic equa¬ 
tion where the negative solution is discarded. The SNR is 
calculated within a 2x2 spaxel extraction aperture, although 
for the 30 x 60 mas scale we upsize to a 120 x 120 mas box. 
We compare R, H and K bands. Eor each band we create an 
input cube containing a point source with a flat spectrum 
(constant Mab). The input cubes are generated with a spa¬ 
tial sampling of 0.4 mas for the 4x4 mas output scale and at 



A [/im] 

Figure 5. Comparison of the total SNR and object signal be¬ 
tween HSIM (red) and the SINFONI ETC (grey). The two main 
panels show the total SNR and object signal within two spaxels. 
The two smaller panels show the respective residuals after median 
smoothing by 10 spectral pixels. 


1 mas for the 10 x 10, 20 x 20 and 30 x 60 mas output scales. 
The input cube is generated with ~ 30 wavelength channels 
so as to be computationally efficient. The SNR is an aver¬ 
age value over all the wavelength channels in the spectrum 
extracted from the 2x2 spaxel box. For each grating we 
choose wavelength regions minimally affected by sky lines 
and telluric absorption. The central wavelengths for each 
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input cube are 0.677/xm for i?-band, 1.76/xm for i7-band, 
and 2.145/xm for X-band. 

These sensitivity calculations are intended to provide 
a general indication of the performance of the current in¬ 
strument design. Whilst the instrument parameters (e.g. 
throughput, dark current, read-out noise) are not set at the 
minimum specification, we use simple 2x2 spaxel aperture 
extraction to present reasonably conservative performance 
estimates. Improved SNR would be obtained using the more 
rigorous method of optimal extraction ( Horne|1986 Robert- 
so n|1986 ). 


7.1 Simulation parameters 


We set the following instrument and site parameters as 
to provide conservative estimates for the resulting perfor¬ 
mance. 

Site and telescope: We use a Paranal-like site (and tele¬ 
scope) temperature of 280.5 K. Seeing is set to the Cerro 
Armazones median value of 0.67'' FWHM at 0.5/xm. We 
use the official 39 m E-ELT model, which has a maximum 
all-glass diameter of 37 m, an obscuration ratio of 0.3 and six 
spider arms, giving a total collecting area of A = 932.5 m^. 
We use an emissivity of 0.244 motivated by the current tele¬ 
scope design of seven warm reflections (5 mirror E-ELT + 2 
further mirrors into HARMONI on the Nasmyth platform) 
of protected silver and aluminium (MgE 2 on Ag+Al) coated 
mirrors and six unmasked spider arms. 

Instrument: We use an instrument throughput of 35 per 
cent. The operating temperature of 120 K means the instru¬ 
ment contributes negligible thermal background. 

Detector model: The NIR detector statistics are modelled 
on the KMOS HAWAII-2RG detectors ( Einger et al.|20()8 ). 
We use a dark current value of 0.0053 e“ s“^ pix“^ and read- 
out noise of 2.845 e“ pix“^. Eor the visible CCD detector 
statistics we adopt a dark current of 0.00042 e“ s“^ pix“^ 
and read-out noise of 2.0 e“ pix“^. The wavelength cut be¬ 
tween visible and NIR detectors is set at 0.8 /xm. 


7.2 Results 

We present the results of sensitivities for LTAO observa¬ 
tions in tables and Table gives the estimated limiting 
magnitudes to achieve SNR = 5 measured from a 2 x 2 
spaxel aperture (4x2 for the 30 x 60 mas scale). These are 
separated into visible R-band, and NIR H and K bands 
and for each resolving power. In the R-band sensitivity in¬ 
creases with spaxel size. AO performance at visible wave¬ 
lengths is poor, with Strehl ratios of ~ 1 per cent, so larger 
spaxels incorporate more object flux. In the R-band the 
20 X 20 mas scale offers slightly improved sensitivity com¬ 
pared to the 10 X 10 mas scale and ~ Imag greater sen¬ 
sitivity compared to the 4 x 4 mas scale. However, in the 
R-band, the 10 x 10 mas is ~ 0.4 mag more sensitive than 
the 20 X 20 mas scale. This is due to the increased thermal 
background in R-band, which penalises the larger spaxel 
scales. We also perform identical simulations for SC AO and 
seeing-limited observations. We find that SC AO results in 
~ 0.6 mag better sensitivity in R-band and ~ 0.4 mag im¬ 
provement in R-band. R-band observations with SCAO are 
not possible due to the visible light being used for AO cor¬ 
rection. Seeing-limited limiting magnitudes are ~ 1.5 mag 


Table 2. Point source sensitivity predictions for HARMONI with 
LTAO calculated from a 2 x 2 spaxel aperture centred on the 
object. 


2x2 spaxel aperture 


Band 

R 

4x4 

10 X 10 

20 X 20 

30 X 60 

R 

3500 

22.7 

23.6 

24.3 

25.1 

R 

3500 

26.2 

27.0 

27.0 

26.4 


7500 

25.2 

26.0 

26.1 

25.7 


20000 

24.1 

25.0 

25.1 

24.9 

R 

3500 

25.7 

26.3 

25.8 

24.8 


7500 

25.2 

25.8 

25.4 

24.5 


20000 

24.2 

25.1 

24.8 

23.9 


Notes: Limiting AB magnitudes to achieve SNR = 5 per spectral 
resolution element for five hours on-source (Texp = 900 s). We 
assume sky subtraction with no penalty to on-source exposure 
time. Band column gives the bandpass region, R is the resolving 
power and the remaining columns give the limiting magnitude for 
each HARMONI spaxel scale. 4x2 spaxel aperture for 30 x 60 mas 
scale. 


Table 3. Point source sensitivity predictions for HARMONI with 
LTAO calculated from a square aperture that encloses 50 per cent 
of the ensquared energy (EE) for each band. 


50 per cent EE square aperture 


Band 

R 

4x4 

10 X 10 

20 X 20 

30 X 60 

R 

3500 

23.6 

24.6 

25.1 

25.4 

R 

3500 

24.5 

25.4 

26.0 

26.2 


7500 

23.5 

24.5 

25.1 

25.6 


20000 

22.5 

23.4 

24.1 

24.7 

R 

3500 

24.7 

25.1 

25.1 

- 


7500 

24.1 

24.6 

24.8 

- 


20000 

23.1 

23.9 

24.2 

- 


Notes: Limiting AB magnitudes to achieve SNR = 5 per spec¬ 
tral resolution element for five hours on-source (Texp = 900 s). 
For the R, R and R bands the 50 per cent EE apertures corre¬ 
spond to 400 X 400, 140 X 140 and 80 x 80 mas respectively. We 
assume sky subtraction with no penalty to on-source exposure 
time. Band column gives the bandpass region, R is the resolving 
power and the remaining columns give the limiting magnitude for 
each HARMONI spaxel scale. 


lower than LTAO for R-band at the 30 x 60 mas scale and 
~ 2.5 mag lower for R and R-bands. 

Table shows magnitudes for identical observations as 
table but after extracting the spectrum from an aperture 
containing 50 per cent of the ensquared energy of the PSF. 
For the R, R and R bands these correspond to 400 x 400, 
140 X 140 and 80 x 80 mas respectively. We omit values for 
the 30 X 60 mas spaxel scale in R-band because the 4x2 
spaxel aperture already contains greater than 50 per cent of 
the ensquared energy. This table shows that a greater SNR 
can be achieved in the R-band by increasing the aperture 
size. For the 4x4 and 10 x 10 mas scales, a factor of ~ 1 mag 
improvement is gained. However, for both R and R bands, 
the limiting magnitudes are lower because the extra spaxels 
from the larger extraction aperture contribute more detector 
noise, as can be seen from Figs and 

To illustrate the different noise regimes. Figs and 
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Figure 6. Heat map showing contributions to the variance for LTAO point source observations in if-band for five hours on source 
(T’exp = 900 s) after extracting in a 2 x 2 spaxel aperture centred on the object. The four quadrants show the contributions to the 
variance from (clockwise from top left) sky, telescope, read-out noise and dark current respectively for each observing mode of the 
instrument (spaxel scale and resolving power). The instrument thermal background is assumed to be negligible. Dark colours show 
a greater contribution and light colours show a smaller contribution. For each instrument configuration, the box with the greatest 
contribution contains its percentage contribution to the total variance. Observations in if-band are predominately read-out noise limited 
except at the coarsest spaxel scale and lowest resolving powers which are sky limited. 
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Figure 7. Same as Fig. I^but for X-band observations: heat map showing contributions to the variance for LTAO point source observations 
in AT-band for five hours on source (Texp = 900 s) after extracting in a 2 x 2 spaxel aperture centred on the object. Comparing to Fig.|^ 
the AT-band observations are dominated by the telescope thermal background, except at finest spaxel scale and highest resolving powers 
which are read-out noise limited. 


show heat maps representing the contributions to the total 
variance from each source: sky, telescope, dark current, and 
read-out noise, when observing with LTAO in the iL-band 
and A"-band respectively. Each observing configuration is 
represented by a single box. Light colours show small contri¬ 
bution and dark colours show larger contribution. For each 
observing configuration, the box with the largest contribu¬ 
tor contains the percentage of the total variance from that 
source. It is clear from these figures that in iL-band HAR- 
MONI is generally detector read-out noise dominated, ex¬ 
cept for the configuration of 30 x 60 mas scale combined with 
the two lowest resolving powers, which are sky dominated. 


However, in the A"-band observations are telescope (ther¬ 
mal) background dominated, except for the finest spaxel 
scale and highest resolving powers. 


8 CASE STUDY: SIMULATIONS OF 

EMISSION-LINE GALAXIES AT Z - 2 - 3 

One of the major science goals of HARMONI is to spa¬ 
tially resolve the interstellar medium (ISM) of high-redshift 
(z ~ 2-5) galaxies and measure the physical processes oc¬ 
curring on scales of individual Hll regions. At these red- 
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shifts, the comoving star-formation density was substan¬ 
tially higher than at z ~0, and so this era has been her¬ 
alded as the peak epoch of galaxy formation when most 
of todays massive galaxies formed the bulk of their stellar 
mass ( Madau fc Dickinson||20l4 ). At these early times, the 
Hubble sequence was not in place, and galaxies appear to 
be undergoing significant changes to their morphologies and 
stellar populations. Over the past decade major advances 
have been made in measuring the dynamics of galaxies at 
this epoch using integral field spectroscopy (e.g. SINS: Gen- 
zel et al.||2006[ ^rster Schreiber et al.||2006 [ |2009 OSIRIS: 


Law et al. |2007J |20Q9| KMOS: [Wisnioski et al. 2015| 

also |Glazebrook|2013| review). However, due to the photon- 
starved nature of the observations, deriving the dynamics 
of high redshift galaxies requires long integration times, and 
only offers limited spatial resolution. Moreover, even with 
natural- or laser-guide star AO assisted observations, spa¬ 
tially resolved studies are limited to > 1 kpc resolution (e.g. 


Forster Schreiber et al.|[2011a|b Swinbank et al.|[2012b|a| ). 

Gravitationally lensed star-forming galaxies do provide the 
unique ability to resolve regions at scales of ~60-200 pc 


(Jones et al. 2010 Livermore et al. 2012 2015), however 


these are notable rare exceptions limited to a handful of 
objects. 

The increased light grasp and spatial resolution of the 
E-ELT should allow the study of high-redshift galaxies on 
the scales of individual star forming Hll regions (down to 100 
pc scales). The properties (e.g. sizes, luminosities, velocity 
dispersions, chemical make-up, and spatial distribution) of 
Hii regions reflect the underlying ISM (such as gas density 
and pressure), which in turn reflect the dominant route by 
which galaxies accrete the bulk of their gas. 

The goals of the simulations presented in this section are 
to investigate the ability to detect and measure kinematics 
of emission-line galaxies at redshifts of z ~ 2-3, using HAR¬ 
MONI on the E-ELT. Our goals are to: 

(i) determine how well the global kinematics (e.g. rotation 
curves) can be derived as a function of star-formation rate, 
size and morphology for galaxies at z ~ 2-3; 

(ii) determine the smallest physical scales for which physical 
properties can be derived, including identifying / measuring 
the properties of individual star forming regions. 

To demonstrate the use of our simulation pipeline in this 
area, we present results from a set of simulated observations 
of a sample of emission-line galaxies at ~2-3. We focus 
on galaxies with prominent Ha emission, which falls into 
the K-band at these redshifts. We highlight that our sim¬ 
ulations are not designed to test galaxy formation models. 
Rather they are designed to test how well measured proper¬ 
ties (e.g. rotation curves, clump properties) can be derived 
for a given exposure time (or at a given spatial resolution) 
using a reasonable input galaxy image/spectrum with a set 
of disc and clump scaling relations. 

In the rest of this section, we discuss how we construct 
mock data-cubes, our simulation runs, analysis methods and 
derived conclusions for this science case. 


8.1 Input data-cubes 

We generate a sample of input galaxy data-cubes covering 
a range of star formation rates, redshifts and morphologies 
ranging from smooth exponential discs of gas and stars to 


‘clumpier’ galaxies (gas follows stars in each case). Each 
galaxy datacube is generated with values randomly picked 
from a range of uniformly distributed physical parameters: 
redshift (2.05-2.8), total disc star formation rate SER (1- 
200 M© yr“^), gas fraction (0.1-0.9), inclination (20-70deg), 
position angle (0-360deg), reddening Ay (0-1.5mag), half- 
light radius rni (0.5-2.5kpc), disc intrinsic velocity disper¬ 
sion (15-40kms“^) and metallicity (Z = 0.05-1 Z©). The 
underlying light profile of the galaxy disc follows an expo¬ 
nential profile, and the velocity field follows a simple arctan 
model (Gourteau 1997). We use the star formation law of 


Kennicutt (1998). We also add a number of star-forming re¬ 


gions using scaling relations inferred from observations of 
lensed star-forming galaxies at z ~ 1-3 (e.g. Jones et al. 
2010 Livermore et al. 2012| 2015). The number of star¬ 


forming clumps is set using the redshift dependent clump lu¬ 
minosity function from Livermore et al. (2012 2015), where 
the normalisation is a function of the disc gas fraction and 
ranges from 0.01-2. The clump velocity dispersion and sizes 
use the scaling relations from Livermore et al.| ( |2015| ) . 

Einally, we include stellar continuum assuming either a 
constant or exponential star formation history (with an in¬ 
tegral that matches the dynamical mass after accounting for 
gas fraction). We use a solar metallicity simple stellar pop¬ 
ulation model with a Ghabrier (2003) initial mass function 
to derive the stellar continuum and assign this to the disc 
according to its luminosity profile. 

Eor the purposes of this analysis, we generate input 
cubes with a spatial sampling of 10 x 10 mas and a resolving 
power of 10,000. 


8.2 Simulation runs 

We simulate a series of mock HARMONI observations of 
these galaxies using LTAO. We adopt Texp = 900 s and use 
the R ~ 3500 H K grating. We focus on simulations at 
20 X 20 mas (~ 200 pc at z = 2-3) scales unless otherwise 
stated. The pipeline parameters (site, telescope, instrument 
and detector) are all set identically to those used for the 
sensitivity calculations (see section [7T| . 


8.3 Global kinematic measurements 

In Eig.j^we show the recovered Ha flux distribution and gas 
kinematics for a sample of five smooth-disc galaxies assum¬ 
ing a 10 hour integration. We use a Gaussian fitting routine 
to fit the Ha (and Nil) emission lines spaxel-by-spaxel. This 
routine iterates over each spaxel and fits both the contin¬ 
uum and Gaussian profiles to the spectrum. In cases where 
no fit is made we average over the surrounding spaxels to 
increase signal at the expense of spatial resolution. Eor the 
velocity maps, which show the global gas kinematics, we bin 
over a 3 X 3 spaxel box, giving a varying resolution of 20 mas 
in the bright regions to 60 mas in fainter regions. We use a 
SNR threshold of 7 for detection of an emission line. 

Erom Eig. I^we see that HARMONI is capable of mea¬ 
suring velocity profiles in galaxies down to Milky Way-like 
star-formation rates on scales of at least ~ 200 pc (in the 
brightest/highest SNR regions) in this integration time. In 
Fig. I^we show the rotation curves for galaxy B (extracted 
along a ~ 1 kpc wide slit aligned along the semi-major axis) 
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Figure 8. Maps of Ho; flux (top row) and line-of-sight velocity (bottom row) for smooth disc galaxies. Galaxies are ordered in increasing 
SFR from left to right. Also shown is the redshift of each galaxy. All galaxies are observed for 10 hours at the 20 x 20 mas (~ 200 pc) 
scale. Velocity gradients are easily measured even for the lowest SFRs with a factor of ~ 5 improvement over existing instruments in the 
highest SNR regions. 



r[kpc] r[kpc] r[kpc] 


Figure 9. Rotation curves for smooth disc galaxy B with varying SFR (increasing left to right) at 20 mas scale (red circles) as observed 
by the E-ELT with HARMONI. The input model curve is shown as the dashed black line. Also plotted in the centre subplot is the 
rotation curve for a 1 hour observation at the 30 x 60 mas scale (blue squares), which has been binned up to 60 mas spaxels. We also show 
a simulated 10 hour VLT (SINFONI) observation at 100 mas (green triangles). All curves are extracted along a ~ 1 kpc slit aligned along 
the semi-major axis. From these simulations the E-ELT with HARMONI offers a 10 fold improvement in observing efficiency compared 
to current telescopes, i.e. the E-ELT could observe 10 times as many objects in equal observing times. It also offers finer resolution, with 
5x more independent data points at the 20 mas scale, and improved sensitivity for equal observing times on the same object 


for the same galaxy properties, but with a star-format ion 
rate of 3.5, 7 and 14 M© yr“^ (increasing left to right). This 
galaxy has the steepest rotation velocity from our sample 
so gives the best indication of which regions can be recov¬ 
ered for a small range of quiescent star-format ion rates. It 
also shows the effect limited spatial resolution can have on 
tracing the inner part of the rotation curve (‘beam smear¬ 
ing’). The velocity profile is traced very closely by HAR¬ 
MONI at 20 mas sampling through the central part of the 
galaxy for all SFRs. The half-light radius of this galaxy is 
rhi = 1.3 kpc so the curve is recovered out to 3rhi for an 
SFR of 3.5 M© yr“^, which increases to 4.6rhi for a four-fold 
increase in SFR. 

In Fig. we also show the rotation curve for a 1 hour 


E-ELT (HARMONI) observation at the coarser 30 x 60 mas 
scale (blue squares), which has been binned up to 60 mas 
spaxels. Comparing this to the rotation curve derived from a 
simulated 10 hour VLT (SINEONI) observation at 100 mas 
we see that the E-ELT offers higher resolution data at better 
sampling with 10 x greater efhciency in observing time. Eor 
equal observing time the VLT curve underestimates the true 
curve at all radii and also only extends to 3rhi, compared 
with 4.2rhi from HARMONI at 20 mas sampling. Deriving a 
simple dynamical mass estimate from each observation gives 
- 9 X 10^ M© from E-ELT (HARMONI) and - 4 x 10^ M© 
from the VLT (SINEONI) simulation. Comparing to the in¬ 
put value of 8.6 x 10^ M© this represents both an accurate 
estimate and a factor of 2 improvement by HARMONI af- 
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ter 10 hours observing. From these simulations we find that 
the E-ELT with HARMONI offers a 10 fold improvement 
in observing efficiency compared to current telescopes, i.e. 
the E-ELT could observe 10 times as many objects in equal 
observing times. It also offers improved sensitivity and finer 
resolution, with 5x more independent data points at the 20 
mas scale, for equal observing times on the same object. 


8.4 Detailed kinematics 

The fine spatial sampling of HARMONI coupled with LTAO 
will allow for very detailed observations of z ~ 2 galaxies. In 
Eig.[^we show the observed Ha intensity and velocity dis¬ 
persion maps for the clumpy galaxies in our simulations. To 
maintain the high spatial resolution required for detecting 
individual star forming regions we bin over a 2 x 2 spaxel 
box where no fit is made to the emission lines, giving an 
effective resolution of 40 mas 300 pc) in fainter regions, 
and again use a SNR threshold of 7 for detection of an emis¬ 
sion line. As Eig.p^ shows, detailed structure is seen in both 
the Ha emission maps and a maps. The maps of galaxy J 

3 spaxels in di¬ 


in Eig. 10 (far right) show structure of 
ameter. Thus HARMONI will be able to make very detailed 
measurements of galaxy substructure, even for galaxies with 
Milky Way star-formation rates. 

To demonstrate HARMONFs ability to discern proper¬ 
ties of individual star forming regions, we focus on the galaxy 
with the lowest star-formation rates from the ‘clumpy’ 


galaxy simulations (galaxy E in Eig. 10). Eig. 11 shows the 
recovered Ha signal-to-noise maps for a 10 hour observa¬ 
tion at the 20 x 20 mas (top row) and 10 x 10 mas (bottom 
row) scales. Eor star-formation rates below lOM0yr“^ the 
20 mas scale offers superior signal-to-noise, while still resolv¬ 
ing individual clumps. However, at a star-formation rate of 
lM0yr“^ no clumps are detected with a SNR > 5 at the 
20 mas scale in 10 hours. 

To quantify the likely detection and properties of 
the clumps in a HARMONI observation we use the 2D 


CLUMPFIND routine by Williams, de Geus & Blitz (1994) 


to determine positions and sizes of individual clumps. This 


method has been used previously by Livermore et al. (2012 


2015) on observations of lensed galaxies. The routine uses 


isophotes to define clumps starting in the brightest regions 
and then moving down through the isophote levels. Any iso¬ 
lated contours are defined as new clumps, and any which 
enclose an existing peak are allocated to that clump. A con¬ 
tour which encloses two or more existing peaks has its pixels 
divided between them using a ‘friends-of-friends’ algorithm. 
We follow a similar procedure as used in [Livermore et al.| 
(2015) and set the minimum threshold at 3a and move up 


in Icr steps. 

In galaxy E we detect three clumps (indicated by yel¬ 
low circles on Eig. 11), and proceed to calculate their lu¬ 


minosities and then star-format ion rates by summing the 
pixels within the clump from the Ha map (corrected for 
local background from the underlying disc). We measure 
star-format ion rates in these three clumps to be 1.46 ±0.02, 
0.26 ± 0.01 and 0.23 ± O.OlM0yr“^, where the uncertain¬ 
ties are derived from the Ha variance map. As a crosscheck 
of our method we integrate the complete observed Ha map 
and calculate the total star-format ion rate of the galaxy to 


be 4.12 ± 0.03 M0 yr“^ (which is similar to the input value 
of 4.35 M0 yr“^). 

Einally, assuming the clumps have circular symmetry, 
we infer a radius of each clump. We measure the radii as 
980 ± 150, 500 ± 200 and 520 ± 250 pc respectively. These 
match closely to the input clump sizes of 1000, 630 and 
500 pc respectively. Comparing these measurements to ob¬ 
servations of lensed galaxies (Eig. 9 of Livermore et al.|2015 ) 
we see that HARMONI will be capable of detecting and 
measuring properties of clumps at least a factor 2 smaller 
than currently possible for normal (unlensed) galaxies at 
z 2 500 pc compared to ~ Ikpc). Thus, it should 

be possible to observe the same galaxy at the 10 mas scale 
for a greater number of hours and measure properties of 
even smaller clumps. In fact, from the signal-to -no ise map 
of galaxy E with SFR = lO.2M0yr“^ in Eig. we see 
there are clumps of three pixels diameter in the 20 mas scale 
with SNR > 5. This corresponds to ~ 250 pc. Our current 
analysis is a first step to showcase the capabilities of the 
simulation pipeline, and we will undertake a more thorough 
analysis of HARMONFs ability to detect and measure star 
forming clumps in a follow-up paper. 


9 CONCLUSIONS 

We present HSIM: a new simulation pipeline for the HAR¬ 
MONI integral field spectrograph on the E-ELT. The 
pipeline takes input data-cubes and simulates observations, 
folding in sky, telescope, instrument and detector param¬ 
eters to create output mock data. HSIM is able to provide 
quantitive measures of the precision with which we can de¬ 
rive a number of key physical parameters for particular sci¬ 
ence cases. It allows the user to gain an understanding of the 
uncertainties associated with making a specific astrophysi- 
cal measurement. We have described HSIM and presented two 
studies: point source sensitivity estimates and simulations of 
z ~ 2-3 emission-line star forming galaxies. Our main con¬ 
clusions are as follows: 

1. The E-ELT AO PSE is predicted to be a strong func¬ 
tion of wavelength. We show the importance of incorporating 
this into our simulations. We use a novel parameterisation 
method to create specific PSEs at each wavelength of the 
data-cube, which then is convolved with its corresponding 
data-cube channel. 

2. HSIM has been thoroughly crosschecked with exist¬ 
ing implementations for other instruments, including the 


pipeline of Puech, Yang & Elores (2010) and the ESO SIN- 


EONI exposure time calculator, and these are consistent 
with our code. 

3. We derive point source sensitivity estimates for ob¬ 
servations with HARMONI using LTAO. The 20 mas scale 
offers the greatest point source sensitivity in iF-band, but 
the 10 mas scale is more sensitive in iF-band due to the in¬ 
creased thermal background. We also show that iF-band ob¬ 
servations are predominately read-out noise limited, whereas 
FF-band is strongly thermal background limited. 

4. We perform simulations of z ^ 2-3 star forming 
emission-line galaxies and find that the E-ELT with HAR¬ 
MONI will be capable of obtaining velocity maps of these 
galaxies down to Milky-Way SERs at these redshifts, with 
a factor ~ 5 improvement in spatial resolution over current 
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Figure 10. Maps of Ha flux, line-of-sight velocity, and velocity dispersion for clumpy galaxies. Galaxies are ordered in increasing SFR 
from left to right. Also shown is the redshift of each galaxy. All galaxies are observed for 10 hours at the 20 x 20 mas scale. Velocity 
structure is evident even at the lowest SFRs and the velocity dispersion maps show internal structure of ~ 5-6 spaxels diameter 500 pc 
radius) for the faintest galaxy down to ~ 3 spaxels 250 pc radius) for the brightest galaxy. 



Figure 11. SNR maps of a model clumpy galaxy for varying SFRs observed at 20 X 20mas (top row) and 10 x 10mas (bottom row). 
SFR increases from left to right and is denoted in the top left of each panel. The contours show constant SNR increasing from 2.5 to 
12.5 in steps of 2.5. The size of the LTAO PSF is shown in the bottom left of the each map in the top row. The total exposure time is 
10 hours. The three star-forming clumps detected by CLUMPFIND are highlighted with yellow dashed circles in the central panel of the 
top row. For comparison we also show the input morphology (10 mas sampling) in the right hand subplot. Sub-kpc sized star forming 
regions are detected with SNR > 5 for SFR ~ 4 Mq at 20mas. Similar regions are detected for SFR ~ 10 Mq at lOmas. 
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generation instruments. By deriving rotation curves we find 
that the E-ELT with HARMONI offers improved resolution 
data with a 10 fold improvement in observing efficiency com¬ 
pared to current telescopes, i.e. in equal observing times the 
E-ELT could observe 10 times as many objects. Eor equal 
observing times on the same galaxy the increased sensitiv¬ 
ity of the E-ELT with HARMONI gives a factor of ~ 2 
improvement in dynamical mass estimates. 

5. HARMONI with LTAO will provide exquisite re¬ 
solved spectroscopy of z ~ 2-3 galaxies, allowing the detec¬ 
tion of individual star forming complexes and measurements 
of their properties. We demonstrate that it will be possible 
to detect star forming Hll regions down to at least ~ 500 pc 
radius in a galaxy of SER = 4 M© yr“^ in a single night. We 
stress that this is only a representative example and a more 
detailed analysis of star forming clumps will be presented in 
a later work. 
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